function [ Fl ] = Flfcn(alfa,s,micro,Fz)
%FL Summary of this function goes here
%   Detailed explanation goes here
a1=-21.3;
a2=1144;
a3=49.6;
a4=226;
a5=0.069;
a6=-0.006;
a7=0.056;
a8=0.486;

D=a1*Fz^2+a2*Fz;
C=1.65;
B=(a3*Fz^2+a4*Fz)/(C*D*exp(a5*Fz));
E=a6*Fz^2+a7*Fz+a8;
psi=(1-E)*s+(E/B)*atan(B*s);

Fl=D*sin(C*atan(B*psi));

end

